UNCLASSIFIED 


AD  NUMBER 


AD837464 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  JUL  1968. 
Other  requests  shall  be  referred  to 
Department  of  the  Army,  Fort  Detrick,  MD 
21701. 


AUTHORITY 


SMUFD  D/A  ltr,  14  Feb  1972 


THIS  PAGE  IS  UNCLASSIFIED 


V 


r 


SO 

CO 


TRANSLATION  no. 


DATE:  S~*  LY 


/rtf 


DDC  AVAILABILITY  NOTICE 

Reproduction  of  this  publication  in  whole  or  in  AUG  1  - 1968 

part  is  prohibited.  However,  DDC  is  authorized  .  - 
to  reproduce  tie  publication  for  United  States  LJiji _ ; 

Government  mtrnnfipfi; 


tt  A.  vovtnjointu 

srrss^/^rir^.-Mr* 

5,l"'c?'/TT»‘  Technical  Rciease  B^ch/ 
ID,  Frederi eft,  Maryland  2! 701 


DEPARTMENT  OF  THE  ARMY 
Fort  Detrick 
Frederick,  Maryland 


Best 

Available 

Copy 


HEATING  AND  COOLING  OF  SIMPLE  GEOMETRIC  BODIES 


by  Dr,  ling,  Heinrich  Groeber,  Berlin-Wilmeradorf 

Communication  from  tlie  Heat  Propagation  Committee  of  the  Association  of 
German  Engineers,  (This  £83ay  will  Portly  be  published  in  an  expanded 
version  aa  researoh  report.) 

mwmwMmmmiiwiMi/ym 

doit.  Vereinos  Doutsch.  Inronieurs  69(21):  709-11,  1??5. 

This  article  features  graphs  and  curves  on  the  heating  and  cooling 
time  of  spherical  bodies,  cylinders,  and  plates  so  that  these  developments 
can  be  easily  traced.  In  looking  at  these  results,  we  must  keep  in  mind 
that  in any  tasks  arising  In  furnace  and  heating  technology,  metal-working 
and  metal  improvement  can  bs  traced  back  to  the  heating  and  cooling  of 
these  simple  bodies.  In  some  cases,  the  technical  literature  on  the  sub¬ 
ject  does  not  offer  any  data  on  the  heat  transfer  /propagation/’  numbers. 
Here  we  can  use  the  graphs  in  order  to  find  the  characteristic  quantities 
and,  fi’om  this,  the  heat  transfer  number,  a  supplementary  note  mentions 
a  work  by  Williamson  and  L,  H,  ,idams. 

The  technical  theory  of  heat  transfer  has  in  the  past  been  concerned 
mostly  with  heat  transfer  /conduction^  processes  in  the  settled  state,  al¬ 
though  the  processes  of  heating  and  cooling  of  solid  bodies  play  a  major 
role  eveiyvhere  in  the  bast  area  of  technology. 

•Mathematical  phyeics  has  been  exploring  these  processes  for  the 
simplest  bodies  quite  accurately  for  a  long  time  but  the  results  do  not. 
show  up  ir.  the  technical  literature  on  the  subject.  The  reason  for  this 
is  to  be  found  in  the  unwieldiness  and  confusion  of  the  final  formulas, 
most  of  which  contain.  Fourier  or  other  similar  endless  series.  But  the 
practical  operator  in  moat  instances  does  not  have  the  necessary  background 
nor  does  he  have  the  time;  he  cannot  calculate  and  compute  for  days  on  end, 
just  to  get  a  single  numerical  value. 

This  being  the  situation,  it  might  be  a  good  idea  to  have  one  man 
tackle  the  computing  and  calculation  job  and  to  work  out  several  important 
problems  for  the  most  frequently  encountered  numerical  ranges.  This  would 


then  give  the  general  public  and  the  specific  user  an  easily  handled  body 
of  numerical  data.  The  management  group  of  the  Committee  of  Technical  and 
Economic  Experts  of  the  Reich  Coal  Council  for  Fuel  Consumption  and  of  the 
Scientific  Advisory  Council  of  the  Association  of  German  Engineers  immediately 
approved  a  suggestion  along  these  lines  and  together  prevailed  upon  their 
boards  to  make  the  nocessary  fbnds  available,  Mathematics  research  stu¬ 
dent  Heinz  Winterfeldt  than  made  all  of  the  confutations  qn  which  I  would 
like  to  report  below. 

Cooling  the  Sphere 

A  sphere  with  a  radius  of  R,  which  was  first  heated  in  its  entire 
mass  to  a  terrperature  of  ti°C,  is  suddenly  placed  in  the  vicinity  of  lower 
temperature  t2°C.  lhe  sphere  will  then  give  off  its  warmth  to  the  surround¬ 
ing  area;  the  outer  layers  will  cool  off  Quickly  and  strongly  while  the 
inner  layers  will  follow  slowly  until,  at  the  end  —  theoretically  after 
an  infinite  period  of  time  —  the  entire  sphere  has  been  cooled  to  to°C. 

This  entire  process  is  so  very  simple  that  anyone  without  the  slightest 
physics  background  can  understand  it:  this  is  vhy  we  might  be  inclined  to 
think  that  it  should  not  be  too  dif.fa.cult  to  follow  the  whole  process 
mathematically.  The  following  pa&es  will  show  that  this  i8  not  the  case. 

When  we  say  that  we  want  to  follow  the  process  mthematically,  we 
mean  that  we  want  the  mathematical  answer  to  the  following  three  questions: 

1.  What  is  the  time  curve  of  the  surface  temperature? 

2.  What  is  the  time  curve  of  the  temperature  at  the  center? 

3.  What  is  the  time  curve  for  the  heat  loss? 

We  can  simplify  the  way  we  write  all  of  these  equations  from  the  very 
outset  by  figuring  the  temperature  not  in  degrees  Centigrade  or  in  degrees 
of  absolute  division  but  rather  by  setting  the  environmental  temperature 
at  zero.  All  temperatures  of  the  sphere  are  then  considered  to  be  over- 
temperatures  and  can  be  ur-dignated  with  £?  ,  (When  W6  heat  a  cold  sphere 
in  a  hot  environment,  we  automatically  get  $  values  with  a  negative  sign.) 

Similarly,  we  are  net  going  to  figure  with  the  heat  content  of  the 
sphere  but  rather  only  witb  the  heat  surplus  above  the  environmental  tempera¬ 
ture. 


Lot  us  use  the  following  symbols  here: 


R  the  radius  of  the  sphere . . . .  in  m, 

t  the  time.... . . . . . .  in  hrs. 
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e,  the  initial  temperature  of  the  sphere  )  in  degrees  of  the  100- 

©•  the  surface  temperature  of  the  sphere  )  part  scale,  measured  as 

temperature  in  the  center  of  the  sphere  )  overtemperatures. 
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a 

cx 
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thermal  conductivity  of  the  sphere  in  kcal  at  time  0, 
thermal  conductivity  of  sphere  material  in  , 
density  of  sphere  material 


specific  heat  of  sphere  material  in 


kg 

m*’ 

_ Kcal 
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temperature  conductivity*^  of  sphere  material  in  ■  JL. _ 

kcal 

m,0O  h'  ' 


heat  transfer  number 


in 


relative  heat  transfer  number  ss. 
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The  computation  for  the  time  curve  of  the  surface  temperature  then 
give  us  the  expression  (*) 

e.  =  et  y  a  •!"*’* -n00*’’*  9  “ 


*  =  ! 


«lo  9k  CO#  *k 


(la) 


(*)  Derivation  of  this  and  subsequent  formulas:  e.g.,  Groeber,  Die  Grund- 
gesetae  der  V.'aermeleitung  und  cl  es  Vfaermeueberganges  (The  Baeic  Laws 
of  Thermal  Conductivity  and  Heat  Transfer),  1921,  pages  44,  51  and  54, 
Julius  Springer  Publishers. 


Here  of  course  1>Jk  would  be  the  system  of  the  infinitely  large 
number  of  roots  of  the  transcendental  aquation 


r  cos  »  =  (1  —  h  R)  sin  r 

Thera  is  only  one  thing  we  can  easily  tell  from  this  equation:  the 
value  Qn  is  proportional  to  the  value  Q.  •  as  for  the  rest,  the  expreseion 
is  a  rather  confused  function  of  the  four  independent,  variables  t,  R,  a, 
and  h.  Anyone  can  plainly  see  that  the  matheiratical  evaluation  of  this 
equation  is  a  rather  timeconsuming  process  and  that,  no  matter  how  good 
our  mathematical  training,  it  is  impossible  to  get  the  course  /curve7  of 
the  function  from  the  fomula.  Besides,  because  of  the  large  number  of 
independent  variables  it  is  impossible  to  show  the  numerical  values  of  the 
function  in  a  single  numerical  table  or  on  a  single  sheet  of  curves  or 
r  graphs.  In  short,  the  result  of  the  mathemati^l  compution  in  the  form  of 

5:  equation  (la)  has  all  of  the  properties  that  make  it  unsuitable  for  further 

i-:  practical  use. 

2  Now,  there  is  one  Vay  out:  we  can  say  that  is  the  only 

i  variable  /&  single,  unique  variable/ j  and  at  the  cams  time  we  would  of 


course  note  that  the^  values,  for  their  part,  depend  only  of  the  single 
variable  hR.  In  this  way  we  can  conceive  of  the  inf inite ' series  as  a 
function  with  only  two  variables  and  equation  (la)  wiH  then  look  like 
this: 

■*•“*•**(«»•  **).  (lb) 


The  values  h  B 

we  now  call  the  characteristic  quantities;  these  are  always  pure  numerical 
values,  i.e.,  magnitudes  without  dimension.  The  possibility  If  grouping 
the  four  variables  into  groups  of  characteristic  quantities  can  be  ex¬ 
plained  with  the  help  of  the  principle  of  similarity  /similitude/  or  the 
theory  of  dimensions. 

The  numerical  evaluation  of  the  function  $ g  in  equation  (lb),  with 
its  two  variables,  still  requires  a  rather  big  Computation  effort;  but 
once  we  have  done  this  job,  we  can  bring  out  the  results  in  a  single  numeri¬ 
cal  table  or  in  a  single  graph. 


In  a  similar  manner  we  can  also  convert  the  equation  for  the  tempera¬ 
ture  of  the  center.  Here  we  have:  -  _ _ _ 
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Surface  Temperature  of  Sphere  When 


0. 

=  1»: 

©0 

=0. 

*o| 

•  t  ' 

■e 

• 

0.000l|  M01 

04)1 

o.t 

'  1  * 

10 

Ll. 

C'' 

0 

J—  0.01 

1/W 

1,00 

0,99 

0,94 

0,89 !  0,64 

0,88 

0,22 

0,10 

0,00 

0,06 

_ 

1.00 

0.<K> 

O.Rfi 

073 

nan 

n  17 

ft  Aft 
V)V» 

ClGo 

G,00 

n  in 

1  lOv) 

Giujj 

0*60 

0,79 

0/14  0,28 

o,iu 

0,06 

0,02 

o/w 

C,28 

— 

1,00 

0.99 

0,03 

0,64 

0,4/  an> 

aoa 

0/)l 

0,00 

— 

,  M 

— 

l/X) 

0,08 

0,86 

0.46 

0.24  0/H 

o/w 

— 

— 

— 

'  l/> 

— 

1,00 

0,97 

0,73 

0,23 

0,07  0/» 

-  — 

— 

— 

—  • 

I/)0 

0,99 

M3 

0r63 

ao3 

0/» 

— 

— 

— 

_ 

— 

6,0 

l/w 

0,99 

0.23 

MO 

— 

■ — 

— 

— 

_ _ 

10/) 

1,00 

0,97 

0,76 

(MW 

O/W 

— 

— 

— 

— 

•  — 

26/) 

IjOO 

0.94 

0^8 

o/w 

— 

— 

— 

— 

—  - 

-  4  - 


Nuu.-irical  Table  ? 

Temperature  in  Center  of  Sphere  Y/hen 
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Numerical  Table  3 
Heat  Lose  of  Sphere  When  e,  =  io-' 
C  =  (w fcflj. 
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Numerical  Table  k 

Surface  Temperature  of  Cylirdor  Whan 
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The  third  question  we  &uked  above  —  tiie  tine  curve  for  the  hoat  I033 
—  can  ba  answered  with  the  following  equation.  The  heat  volume,  which  the 
sphere  has  given  off  until  time  t,  is 


-  *  Ti»  .  r  ,  *  Vfi  J  ft1"  ■'t-*'*  <=»»%)*  ( ,  '  ’**  A>) 

®  J  *  rt»  rft  —  *Id  COI  \  ..  ' 


Since  $cis  the  overtemperature  of  t  he  sphere  over  /above7  the  en¬ 
vironment,  we  have 


which  is  equal  to  the  heat  surplus  of  the  sphere  at  time  t  =  0. 

In  keeping  with  what  we  said  earlier,  we  can  now  write: 

'.3b) 
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Numerical  tables  1,  2,  and  3  an]  the  curves  in  figured  1,  ?,  and  3 
illustrate  the  values  of  the  throe  functions  for  tt»  sphere.  Here  the 
abscissa  in  tho  illustrations  is  plotted  on  tho  logaritlimic  scale  so  ttiao 
the  entire  range  of  variables  h  ft  from  0.0001  to  50  can  be  ohown. 


Cooling  tho  Cylinder 

A  very  long  cylinder  with  a  radius  of  ft  is  now  suddenly  brought  non 
a  high  tempo nature  to  a  low  environmental  temperature.  We  want  to  compute 
the  time  curve  of  the  surfuca  temperature  ^  ,  the  temperature  in  the 
axis,  and  the  heat  I033  for  a  section  with  a  length  L  of  thin  cylinder. 


Disregarding  the  cooling  of  the  terminal  surfaces  becauso  of  the 
great  length  of  the  cylinder,  we  have  a  computation  hero  which  is  entirely 
similar  to  that  of  tho  sphere.  The  fact  is  that  we  have  tho  following: 


k  =  w 


i _ 


(Aa) 


where  represents  the  3ystem  of  infinitely  many  roots 
^  J-^  (yY )  hIhTQ  ( Jjl )  and  Jq  as  v/oll  as 


tal  equation 
functions. 
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h  -  w 


e* -e/VYl.  *>*)  S'-'* 
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of  the  transcondon- 
are  two  Bouuol 

(5a) 


-  8  - 


Figure  5.  Tempera ture  of  axis  of  cylinder  when 


iron  fomt.l  «*"  «' 

“  v=(w«)*'P  »«)• 

Figure  6.  Heat  loss  of  cylinder  when  ’  - 


r0‘l .1 


The  numerical  evaluation  of  the  three  equations 

=  •  •  •  •  ' 

em  =  ...  ( 

'  g=(»K0)„ ■/'(*«.“£)  •  •  (, 

is  shown  in  numerical  tables  4i  .*>•  and  6  and  in  figures  4,  5>  and  6 


Unilateral  Cooling  of  Plate 


Let  us  take  a  very  large  plat©  with  thickness  X  and  let  us  insular 
it  completely  on  one  side  while  the  other  three  sides  are  free.  If  this 
plate  is  suddenly  brought  from  a  high  temperature  to  a  colder  environment, 
it  will  gradually  give  off  its  entire  heat  content  to  the  environment 
through  the  noninsula ted  surface.  Now  we  are  going  to  have  to  determine 
the  time  curve  of  temperature  Q §  of  the  free  surface,  the  temperature 
of  the  insulated  surface,  and  the  heat  loss  Q  for  a  section  F  of  the  plate  I 


The  computation  gives  us  the  following 


b  JT» 


COli  6,\ 


h  =  i 


vdier eli/  represents  the  infinitely  many  roots  of  the  equation 


>=*•2 
* = i 
k  ^  to 


2 _ «in  »k 

i*  4-  «ln  ik  co*  lit 


Q-fxcy 9. Vs L  _  ~ 'fc' 

dwJ  l5jk*  -f-  dft  Bin  dh  COS  6k  »  i 


<?a) 


(8a) 

(9a) 


The  numerical  values  from  these  three  equations 

e.  =  ©.  .......  .  ^ 

em^e.9m'[h  at,  . ).  (8b) 

..../•  (9b) 


are  shown  in  numerical  tables  7,  8,  and  9  arid  in  figures  7,  8,  and  9. 
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value  a  for  the  functions  if).  £  . 
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The  time  is  put  in  terms  of  hours: 
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Numerical  TabJe  7 

Temperature  of  Free  Plate  Surface  .hen 
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numerical  Table  8 

Temperature  of  Insulated  Plate  Surface  V/hen 
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Numerical  Table  9 

Heat  Loss  of  Unilaterally  Cooling  Plate  V/hen 
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From  the  physical  numerical  tables  we  get  the  foilwing  figures  for 

steel: 


i  =  50  [kcal/m  h  °C], 
y  =  7700  [ke/ms] , 
c  —  0,13  [kc»l/kg»C]; 


from  which  we  can  calculate : 

o«=  — —  60 

e  y  O.lTTWo  =  0,06  I™’/11]  ■ 
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Figure  11.  Cooling  speed  for  center  or  axis  of  various  bodies 

according  to  Williamson  and  Adams.  Legend:  a  -  plate; 
b  -  square  beam,  length  infinite;  o  -  cylinder,  length 
infinite;  d  -  sphere;  e  -  cylinder,  length  equals 
diameter;  f  -  cube. 

Computation  of  characteristic  quantities  from  these  values: 


ft«=- 


•*“?!?• 0.1  =1,0, 


60 


( «>) ,  ~  °oM  ’ 0,01  “  0,06 1 
l** 

kl\  =  0.05.02 

1«Vj  0,01  0,2 


u*), 

on* 

•0,20=1,00. 


I 

I 

I 


V.re  can  now  read  off  directly  and  get  the  following: 


Numerical  Table  10 


(cljWorto  mu»  Zwock* 

,  dor  Uorechnun* 

<i  =  96« 

<•  —  0  mi  a 

<•=  Umin 
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rt 

v»  * w 

n  a 

0,07 

:  2.  der  Tcmporatur 
dor  Mitte . 

_ 

e;_ 

1,00 

0,69 

0,11 

3.  der  obgogebonon 
WSrmo . 

Q 

0,12 

0,47 

0,92 

4.  dor  zurUckgeblie- 
benen  Wkrme  .  . 

i— « 

0,83 

033 

036 

Legend:  a  -  values  for  computation  purposes;  1.  of  surface  temperature; 
2.  of  temperature  in  the  middle;  3.  of  heat  given  off; 

4.  of  remaining  heat. 
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Figure  10.  How  to  Illustrate  the  Determination  of  the  Temperatures 
Along  the  Radius  of  a  Sphere. 

Plotting  the  temperature  distribution.  Vie  plot  the  radius  as  the 
abscissa;  as  the  ordinates,  we  plot,  above  the  terminal  points  of  the 
radius,  on  the  one  hand,  the  three  values  Oc  !  &<-  anc*>  on  °^er 
hand,  the  three  values  Q^/  Ec.  .  These  six  Values  alone  give  us  a 
pretty  good  idea  of  the  attenuation  of  the  temperatures. 

By  drawing  the  illustration  in  Figure  10  we  can  get-  more  detailed 
information  on  the  temperature  curve  along  the  radius.  In  the  center,  this 
curve  must  have  a  horizontal  tangent  for  symmetry  reasons.  A  ong  the 
surface,  however,  the  tangents  of  all  curves  must  go  through  the  same  point 
which  lies  on  tic  abscissa  axis  at  a  distance  of  1/h  JmJ  outride  the  sur¬ 
face  (this  applies  equally  to  ths  cylinder  add  the  plate;  for  explanation, 
see  Groeber,  loc.  cit.,  page  It).  With  the  help  of  these  tangents,  we  can 
then  draw  the  temperature  curves  with  the  kind  of  accuracy  that  will 
suffice  for  technical  Engineering  purposes. 

Computation  of  heat  loss  of  sphere.  The  heat  surplus  of  the  sphere 
l [cooled 7  from  280°  C  via  the  oil  bath  temperature  of  30“  C  is 

=  '  0,1 3  •  7700  •  250  =  104&  [kcal  J .  ! 

o  lUvv 
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Of  this  amount,  the  fraction  given  in  line  3  of  Numerical  Table  10  was 
transferred  to  the  oil  at  time3  to,  to,  and  to;  tlie  rest  remained  in  the 
sphere. 

Supplementary  Comment 

While  thi3  article  was  at  the  printer's.  Prof.  K.  Jakob,  Berlin, 
told  me  about  an  american  work  by  Williamson  and  Adams  who  were  after  a 
similar  objective  (cf.,  Williamson  and  L.  H.  /idam3,  "Temperature  Distri¬ 
bution  in  Solids  during  Heating  or  Cooling,"  Physical  Review.  Vol  XIV 
Series  II,  1919,  page  99).  Unfortunately,  the  authors  always  assume  that 
it  is  not  the  environmental  temperature  but  directly  the  surface  tempera¬ 
ture  which  is  given,  respectively,  changed.  This  is  the  same  thing  as  if 
in  this  work  here,  only  the  perpendicular  series  were  given  for  h  R  and 
h  X-oO.  In  the  first  part  of  the  work  the  authors  furthermore  assume 
that,  after  the  settled  state  has  set  in,  the  surface  temperature  rises 
constantly  from  time  t  =  0  on  in  accordance  with  the  law  Qc  —  const,  t. 
Since  this  case  is  of  no  consequence  from  the  technical  ^engineering 
viewpoint,  wc  need  not  discuss  it  any  further  here. 

In  the  second  part  the  authors  assume  that,  after  the  settled  state 
has  set  in,  the  surface  temperature  suddenly  jumps  by  t-ho  value  and 
that  tho  heating  or  cooling  process  new  takes  placo.  They  accurately 
computed  the  example  of  the  sphere  and  give  the  values  &  /  for 

varying  distance  r  from  the  center  of  the  sphere;  seo  Numerical  table  11. 
The  values  for  r/R  -  0  of  course  correspond  to  the  values  /  O  ^  for 

h  Rs=  oo  in  Numerical  Table  2.  On  the  other  hand,  the  values  of  tho  other 
perpendicular  series  constitute  an  important  supplement  here. 

But  I  think  that  Numerical  Table  12  is  particularly  noteworthy  here; 
this  table  is  likewise  taken  from  that  work.  It  seems  that  the  authors 
have  computed  and  compiled  the  cooling  speed  of  the  center  and  the  axis 
for  the  various  bodies  listed  at  the  top  of  the  table.  See  Figure  11. 

•  Numerical  Table  11 

Temperature  in  the  Interior  of  tho  Sphere  ’When 
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Numerical  Table  12 

Temperatures  at  Center  or  axis  of  Various  Bodies  When 
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Legend:  a  -  plate ;  b  -  sqm re  beau,  length  infinite;  c  -  cube; 
d  -  cylinder,  length  infinite;  e  -  cylinder,  length 
equals  diameter;  f  -  sphare. 

a  comparison  between  the  sphere  and  the  cube  or  between  the  cylinder 
of  infinite  length  and  the  cylinder  whose  length  equals  the  diameter  would, 
I  think,  be  rather  interesting  here.  Although  these  relations  apply  only 
to  an  infinitely  large  heat  transfer  number,  i.e.,_only  for  h  R  —  oo, 
they  nevertheless  do  permit  a  transfor  ^/conversion/  to  processes  with  a 
finite  heat  transfer  number,  if  we  are  careful  enough. 
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